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A multiple scattering formulation is used to calculate the force, arising from fluctuating scalar 
fields, between distinct bodies described by (5-function potentials, so-called semitransparent bodies. 
(In the limit of strong coupling, a semitransparent boundary becomes a Dirichlet one.) We obtain 
expressions for the Casimir energies between disjoint parallel semitransparent cylinders and between 
disjoint semitransparent spheres. In the limit of weak coupling, we derive power series expansions for 
the energy, which can be exactly summed, so that explicit, very simple, closed-form expressions are 
obtained in both cases. The proximity force theorem holds when the objects are almost touching, 
but is subject to large corrections as the bodies are moved further apart. 
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Multiple scattering methods for calculating Casimir 
(quantum vacuum) energies between bodies date back 
to the famous papers of Balian and Duplantier [H, 0, Q . 
More recently, Emig and collaborators [1, Q have pub- 
lished a series of papers, using closely related methods, 
to calculate numerically forces between distinct bodies, 
starting from periodically deformed ones Q. The meth- 
ods were developed independently in papers by Bulgac, 
Marierski, and Wirzba 0, S, who, using the Krein 
formula [10|, obtained results for the interaction of two 
spheres, or a sphere and a plate (for Dirichlet bound- 
ary conditions), and by Bordag who rederived the 
methods developed in Refs. 0,01 in a modified form and 
applied them to calculate the first correction beyond the 
proximity force approximation (PFA) for a cylinder near 
a plane. In Ref. |1^] Bordag generalized the method of 
Ref. Ill to the case of a semitransparent cylinder next to 
a plane. Dalvit et al. 13|, used the argument princi- 
ple to calculate the interaction between conducting cylin- 
ders with parallel axes when one cylinder is inside the 
other. Recently, there appeared papers concerning "ex- 
act" methods of calculating Casimir energies or forces be- 
tween arbitrary distinct bodies by Emig, Graham, Jaffe, 



and Kardar [15|, ll6|. See also Refs. [17|, ll8|. Most ex- 
plicitly, an earlier drafted paper by Kenneth and Klich 
[19| appeared which shows that the basis of the latter 
approach lies in the Lippmann-Schwinger formulation of 
scattering theory [20| . 



We might also mention the worldline numerical method 
of Gies and Klingmiiller (2l|, 0, [H, 111] , but that method 
lies rather outside our discussion here. The same ap- 
plies to the work of Capasso et al. [25|, who calculate 
forces from stress tensors using the familiar construction 
of the stress tensor in terms of Creen's dyadics 26l 27|. 
using a numerical engineering method: Finite-difference 
frequency-domain methods are employed in two dimen- 
sions to obtain forces between metal squares and plates 
to 3% accuracy. 



We will now proceed to restate the multiple scattering 
technique, in a simple, straightforward way, and apply 
it to various situations, all characterized by (5- function 
potentials. 

The general formula for the Casimir energy (for sim- 
plicity here we restrict attention to a massless scalar field) 
is [2| 



E = — Tr In G ^ — Tr In GGn 



2r 



2t 



(1) 



where r is the "infinite" time that the configuration ex- 
ists, and G is the Green's function in the presence of a 
potential V satisfying (matrix notation) 



{-d^ + V)G ^1, 



(2) 



subject to some boundary conditions at infinity. (Details 
will be supplied elsewhere |23].) In the second form of 
Eq. ([1]) we have subtracted the energy of the vacuum, 
by inserting the free Green's function Gq, which satis- 
fies, with the same boundary conditions as G, the free 
equation 



S'Go 



(3) 



Now we define the T-matrix (note that our definition of 
T differs by a factor of 2 from that in Ref. [l5| ) 



T = S -l = V{l + GoVy 



(4) 



We then follow standard scattering theory [20|, as re- 
viewed in Kenneth and Klich 'l^. If the potential has 
two disjoint parts, V = Vi + V2, it is easy to derive the 
following general expression for the interaction between 
two bodies (potentials): 



E12 = -— Trln(l - G0T1G0T2), (5) 



where 



T,^V,il + GoV,)-\ 1 = 1,2. 



(6) 



2 



This form is exactly th at g iven by Emig et al. [l5|, and 
by Kenneth and Khch [l9|. In passing from Eq. ([T]) to 
Eq. (O we have removed self-action terms that are gen- 
erally divergent, but that do not refer to the separation 
between the bodies. 

Elsewhere [2^ we will show that this formulation al- 
lows us to rederive the Casimir interaction between two 
semitransparent plates, and the self-energy of the semi- 
transparent sphere. 



2 + 1 SPATIAL DIMENSIONS 

We now proceed to apply this method to the interac- 
tion between bodies, starting with a 2 -I- 1 dimensional 
version, which allows us to describe, for example, cylin- 
ders with parallel axes. Let the distance between the 
centers of the bodies be R. Then we perform a Fourier 
analysis of the reduced Green's function, defined by 

Go(R + r'-r)= j ^e^^-^(^-^')<?o(R± + rl - r^), 

(7) 

where the reduced Green's function has the expansion 
(as long as the two potentials do not overlap) 

50 = E ^..('*0e""^/™'(Atr')e-™>'5^,„,(«i?), (8) 



where lo = iC, and — k1 + C,'^ . The Fourier-Bessel 
transform of the reduced Green's function is 



27r 



-Krn-m'inR). 



(9) 



Thus we can derive from Eq. ([5]) an expression for the in- 
teraction energy (per unit length L) between two bodies, 
in terms of discrete matrices, 

^ = ^ / In det (1 - fhf^t2) , (10) 

where T denotes transpose, and where the t matrix ele- 
ments are given by 



(dr^) / (dr'^)/™(«:r)e-™*/„,,(Atr')e™'*'r. 



Consider, as an example, two parallel semitransparent 
cylinders, of radii a and 6, respectively, lying outside each 
other, described by the potentials 

Vx = \i5{r ~ a), V2 = A2<5(r' - 5), (12) 

with the separation R between the axes satisfying R > 
a + b. It is easy to work out the scattering matrix: 



{tl)mm' — 27rAia(5„ 



/^(kc) 



1 -|- XiaIm{na)Km{Ka) 



(13) 
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FIG. 1: Plotted is the ratio of the exact interaction energy 
(|19p of two weakly-coupled cylinders to the proximity force 
approximation ([20} as a function of the cylinder radius a for 
a = b. 



Then the Casimir energy E is 
E 1 f°° 

— = — / dnntT\n{l ~ A), 

where A = B{a)B{b), in terms of the matrices 

Al a/^,(Ka) 



Br, 



' (a) = Km+ni' (kR) 



1 + \iaIm'{>ia)Kmi{na) 



(14) 



(15) 



As a check, it is easy to reproduce the result derived 
by Bordag [11| for a cylinder in front of a plane, using an 
evident image method. 

In weak coupling, the formula (|14p for the interaction 
energy between the cylinders is 



AiA2a5 
47ri?2 



E 

T 



X / dxxK'^+,^,ix)If,,ixa/R)Ii,ixb/R).{16) 
Jo 

It is straightforward to develop a power series in a/R 
for the interaction between semitransparent cylinders. 
One merely exploits the small argument expansion for 
the modified Bessel functions Im{xa/R) and Im'{xb/R). 
The result is amazingly simple: 



E _ AiaA2& 1 
L ^ 47ri?2 2 ^ \R 

n=0 



(- 

\R 



2n 



where in terms of the binomial coefficients 

2 / T \ 2k 

I n \ I — ^ I n \ 

Pn 



(17) 



(18) 



Remarkably, it is possible to perform the sums [30], so 
we obtain the following closed form for the interaction 
between two weakly-coupled cylinders: 

\n-i/2 

in. — h\ 

1- I 1 III- 



E 

T 



\1a\2b 




R 



(19) 



3 





1.0 




U.B 




0.6 


Si 






0.4 




0.2 




0.0 



0.000 0.002 0.004 0.006 0.008 0.010 
a/R 

FIG. 2: Plotted is the ratio of the exact interaction energy 
(|19l) of two weakly-coupled cylinders to the proximity force 
approximation (|20p as a function of the cylinder radius a for 
h/a = 99. 



We note that in the limit R — a — b = d^O,d being the 
distance between the closest points on the two cylinders, 
we recover the proximity force theorem in this case, 



V{d) 



A1A2 [2ab 1 



d <C a, 5. 



(20) 



327r V R di/2' 
In Figs. [IH2] we compare the exact energy (|19p with 
I 



the (ambiguously defined) proximity force approximation 
(|20)) . Evidently, the former approach the latter when the 
sum of the radii a + b oi the cylinders approaches the 
distance R between their centers. The rate of approach 
is linear (with slope 3/2) for the equal radius case, but 
with slope b'^/Aa^ when a <^ b. 



3-DIMENSIONAL FORMALISM 

The three-dimensional formalism is very similar. 
Again, details will be supplied in Ref . [29*1 . Let us proceed 
to write down the expression following from Eq. ([5]) for 
the interaction between two semitransparent cylinders: 



1 f°° 

^=^j^ dCtrln(l-A), 



(21) 



where the matrix 



Y,Bu"m{a)Bini,^{b) (22) 



is given in terms of the quantities [the three-j symbols 
(Wigner coefficients) here vanish unless l + l' + 1" is even] 



Bwm{a) 



i-'+'V(2Z + l)(2Z' + 



^ lo ; i^m -m ; l + Aia/p+i/2(Ca)^r+i/2(Ca) 



(23) 



For strong coupling, this result reduces to that found by Bulgac et al. [7] for Dirichlet spheres, and recently generalized 
by Emig et al. for Robin boundary conditions. 

For weak coupling, a major simplification results because ot the orthogonality property (I < I'), 



m —TO 



TO —771 



m— — / 

Then the formula for the energy of interaction between the two spheres is 

2 



21" + 1 



(24) 



E = / — ^(2; + l){2l' + l){2l" + 1) 

'0 Wi" 



4i? 




I 



Kl,^y^{x)ll„^{xa/R)ll^y^{xb/R). (25) 



There is no infrared divergence because for small x the 
product of Bessel bmctions goes like 2;^'^'+' ~' and 
I" < I + V because of the triangle property of the 3j- 
symbols. 



Again, it is straightforward to carry out a power series 
expansion in a/R, which turns out to have a simple form 



E = 



\1a\2b ab 
8R~^ ^ n 
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1 " 

— T 

+ 1 ^„ 



2n- 
2m 



(- 

\R 



2n 



2m 



(26) 
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FIG. 3: Plotted is the ratio of the exact interaction energy 
(|27l) of two weakly-coupled spheres to the proximity force ap- 
proximation (|28|) as a function of the sphere radius a for a = b. 
Shown also by a dashed line is the power series expansion (I26|l , 
truncated at n = 100, indicating that it is necessary to include 
very high powers to capture the proximity force limit. 



Once more, it can be recognized as the following closed 
form: 



XiaX2b 

h = in 

I&ttR 



1 



B.) 



1 _ _ 

. V_R B.) ' 



(27) 



Again, when d = R — a 
theorem is reproduced: 



b <^ a,b, the proximity force 



Vid) 



XiX2ab 
WttR 



\n{d/R), d<a,6. 



(28) 



However, as Figs. [3l|4] demonstrate, the approach is not 
very smooth, even for equal-sized spheres. The ratio of 
the energy to the PFA is 



E 
V 



ln[(l + a)V2a] 
In d/R 



d a,b, 



(29) 



for b/a = a. Truncating the power series (pS)) at n = 100 
would only begin to show the approach to the PFA limit. 
The error in using the PFA formula between spheres ([28)1 
outside its range of validity can be very substantial. 



CONCLUSION 

We have used standard multiple scattering techniques 
to calculate the Casimir interaction between two semi- 
transparent (5-function) spheres and between two semi- 
transparent parallel cylinders. When the coupling con- 
stant is weak, we are able to sum the power series ex- 
pansion in a/ R exactly, and obtain a closed form for the 
Casimir interaction energy. This energy reduces to the 
proximity force limit when the bodies are very close to- 
gether. Formulas which extrapolate away from the lim- 
iting proximity force theorem are commonly referred to 
as proximity force approximations; these are ambiguous 
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FIG. 4: Plotted is the ratio of the exact interaction energy 
(|27p of two weakly-coupled spheres to the proximity force 
approximation ([28} as a function of the sphere radius a for 
b/a = 49. 



and in general, the PFAs do a poor job in describing the 
interaction. These exact results represent the first known 
exact closed-form results for the Casimir interaction be- 
tween two bodies which are not plane surfaces. More 
details and examples will be given in Ref. [2^. 
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